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ABSTRACT 


The  feasibility  of  simulating  a  simple  fLorentz  plasma  with  an 
artificial  dielectric  constructed  from  a  random  distribution  of  in¬ 
ductively  loaded  dipoles  is  investigated.  Both  the  differences  and 
similarities  between  the  macroscopic  electromagnetic  properties  of  the 
artificial  dielectric  aha  a  Lorentz  plasma  are  discussed. 
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mental  data  discussed  in  this  report.  We  would  also  like  to  thank 
Mr.  John  Dietz  for  his  many  helpful  suggestions. 
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1.  INTRODUCTION 


The  complete  simulation  of  the  electromagnetic  pulse,  EMP,  associated 
with  a  nuclear  burst  involves  a  number  of  interrelated  parameters.  The 
most  obvious  parameters  are  the  waveform  and  the  magnitude  of  the  EMP. 
However,  if  the  system  being  subjected  to  the  nuclear  EMP  is  in  a  pre¬ 
ionized  region  such  as  the  ionosphere,  we  would  also  have  to  simulate 
the  ionozation  or  dielectric  constant  of  the  medium  surrounding  the  sys¬ 
tem  such  as  a  missile  or  satellite.  In  this  work  we  investigate  the 
feasibility  of  using  an  artificial  dielectric  consisting  of  a  random 
distribution  of  inductively  loaded  short  dipoles  to  simulate  the  macro¬ 
scopic  electromagnetic  properties  of  a  simple  Lorentz  plasma.1*2 

During  the  past  twenty-five  years  a  number  of  so-called  artificial 
dielectrics  consisting  of  regularly  spaced  rods,  parallel  plates,  metal 
spheres,  etc.  have  been  devised  to  reproduce  the  essential  macroscopic 
properties  of  a  dielectric.  The  ordinary  artificial  dielectric  consists 
of  discrete  metallic  or  dielectric  particles  or  lattices  of  macroscopic 
size.3-7  These  artificial  dielectrics  were  first  actually  conceived  as 
large-scale  macroscopic  models  of  microscopic  crystal  lattices.  The 
practical  motivation  for  the  development  of  the  first  artificial  die¬ 
lectrics  was  the  desire  to  obtain  relatively  inexpensive  lightweight 
materials  that  could  be  used  for  microwave  and  radar  lenses.  Several  of 
the  artificial  dielectrics  proposed  for  microwave  lenses  have  the  macro¬ 
scopic  electromagnetic  properties  of  a  plasma.2*4 » 8» 9> i0 


A  review  of  the  literature  reveals  a  number  of  papers  on  the  plasma 
simulation  properties  of  artificial  dielectrics  consisting  of  a  rigid 
cubic  lattice  of  three-dimensional  gr?.ds . 2 > 4 » 8  Experimental  work,  re¬ 
ported  in  the  literature,  has  been  done  only  on  two-dimensional  grid, 
lattice  structures.2* 8,1 4  That  a  periodic  grid  structure  would  produce 
band  structure  resonances  analogous  to  Bragg  scattering  has  been  dis¬ 
cussed  from  a  theoretical  viewpoint,  but  has  not  been  experimentally 
investigated. 11 » 12 » 13» 


The  rigid  cubic  grid  lattice  structure  has  the  disadvantage  of 
being  somewhat  unwieldy  and  difficult  to  support  and  fit  around  an  ob¬ 
ject  with  curved  surfaces.  In  general,  the  cubic  lattice  grid  structure 
lends  itself  conveniently  to  simulation  of  flat,  slab-like  plasma 
sheaths. 


In  this  paper  we  describe  another  approach  to  the  simulation  of  a 
Lorentzian  plasma.  We  propose  a  granular  pellet-like  artificial  die¬ 
lectric  consisting  of  a  random  distribution  of  styrofoam  spheres  con¬ 
taining  inductively  loaded  dipoles  which  can  yield  macroscopic  electro¬ 
magnetic  constitutive  relationships  similar  to  those  of  a  plasma. 

A  Lorentz  plasma  can  be  represented  from  the  viewpoint  of  macro¬ 
scopic  electromagnetic  theory,  as  a  lossy  dielectric  with  a  complex 
dielectric  constant  given  by1*2 
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e  =  eo  fl  -  m2/(v2  +  to2)  +  j<o2(v/io)/(v2  +  to2)],  (1) 

P  L  p  p  J 

where  eo  is  the  dielectric  constant  of  free  space,  to  is  the  frequency  of 
the  electromagnetic  field,  top  is  the  plasma  frequency,  and  v  is  the 
collision  frequency  between  electrons  and  the  gas  molecules.  We  propose 
to  simulate  this  macroscopic  dielectric  constant  with  an  artificial 
dielectric  consisting  of  a  random  distribution  of  inductively  loaded 
dipoles  encased  in  styrofoam  pellets.  A  schematic  drawing  of  such  a 
dipole  is  shown  in  figure  1.  The  quantities  22,  L,  and  R  are  the  length, 
inductance,  and  resistance  of  the  loaded  dipole  respectively.  As  we 
will  shew,  using  simple  quasi-static  arguments,  if  the  effective  in¬ 
ductance  coL  is  much  greater  than  the  effective  capacitance  of  the  dipole, 
the  permittivity  of  a  random  distribution  of  inductively  loaded  dipoles 
is 


e 

P 


4.  N 
3  eo  L 


<!>' 


+  u* 


+  j 


,R.  1 
3  eo  L  V  m 


<£>  +o,2 


(2) 


where  N  is  the  number  density  of  the  dipole  pellets.  One  can  see  by 
comparing  equation  (1)  with  equation  (2),  that  we  can  simulate  a  plasma 
with  plasma  frequency  ojp  and  collision  frequency  v  by  setting 


4  2£  N 
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Figure  1.  Inductively  loaded  dipole. 


2.  THEORETICAL  DISCUSSION 


*\ 


2.1  Quasi-Static  Approach 

We  will  first  approach  the  problem  of  our  proposed  artificial  die¬ 
lectric  from  a  quasi-static  viewpoint  assuming  that  the  inductive  load 
of  the  dipole  is  much  larger  than  the  driving  point  capacitive  impedance 
of  a  short  dipole.  In  a  subsequent  section  we  will  explicitly  consider 
the  dipole  capacitance.  First,  we  discuss  the  simulation  of  a  tenuous 
plasma  in  which  the  collision  frequency  between  free  electrons  and  mole¬ 
cules,  v,  is  negligible. 

We  want  to  construct  an  artificial  dielectric  from  a  pellet-like 
medium  that  has  an  index  of  refraction  given  by2 

n  =  (1  -  w2/io2)  l/2  (5) 


and  an  intrinsic  impedance  given  by 


where  Sp  has  a  frequency  dependence  of  the  form1*2 

ep  =  e o  ( 1  -  w2/w2)  (7) 


and  where  (op  is  the  plasma  frequency. 

In  general  the  index  of  refraction  of  an  artificial  dielectric  may 
be  calculated  in  a  manner  analogous  to  that  employed  in  calculating  the 
index  of  refraction  of  a  molecular  medium.  Assuming  that  the  random 
obstacles  are  not  too  closely  packed  so  that  we  do  not  have  to  resort  to 
the  Clauslus-Mossotti  relation,  the  index  of  refraction  is  given  by5 

n  =  fsx  (1  +  Nx 


,/e0>  Pj  (1  +  Nxm/p0)j  1 


/  2 


(8) 


where 

N  =  Number  of  scattering  obstacles  per  unit  volume, 

s 

Xe  =  Electric  polarizability  of  a  scattering  obstacle, 

Xm  =  Magnetic  polarizability  of  a  scattering  obstacle, 

ex  =  1  (for  styrofoam), 

Pi  =  1  (for  styrofoam), 
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eo  =  permittivity  of  free  space, 

PO  =  permeability  of  free  space. 

From  equations  (5)  and  (8)  we  see  that  in  order  to  achieve  our  goal 
we  need  an  artificial  dielectric  constructed  from  obstacles  that  have  a 
negative  electric  polarizability,  Xe*  If  the  artificial  dielectric  is  to 
simulate  a  plasma  over  a  spectrum  of  frequencies,  the  value  of  Xe  oust 
also  be  inversely  proportional  to  the  frequency  squared.  We  would  also 
like  the  vajJue  of  the  magnetic  polarizability  of  the  embedded  obstacles, 
Xm,  to  be  .«  all  or,  if  possible,  zero. 

A  promising  candidate  for  a  scattering  obstacle  is  an  inductively 
loaded  electric  dipole.  Consider  the  very  idealized  inductively  loaded 
dipole  schematically  shown  in  figure  2.  Let  an  electric  field  E  be 
applied  parallel  to  the  line  between  two  spheres  separated  by  the 
distance  21.  The  field  E  then  establishes  a  voltage  V=2£E  e3ut  between 
the  spheres.  Assume  that  inductance  between  the  two  spheres  is  L,  and 
also  assume  that  Loj  >>  1/wC  where  C  is  the  capacitance  between  the  two 
spheres.  Since  we  are  dealing  with  a  time  harmonic  field  assume  that 
the  charge  accumulated  on  the  dipole  spheres  is  given  by 


q  -  Qe 


then 


24E  eja)t  *  L  — 


2q 

dt2 


=  LQ(-w2)e^Wt 


such  that 

Q  ~  -  2£E/Lw2 
or 

q  =  -  2£E  e^Ut/Lw2. 

The  dipole  moment  of  the  inductive  dipole  in  figure  2  is  then 
p  =  2Zq  «  -4£2E  ej<0t/Lo)2 


or 


XeE  e 


jut 


(9) 


(10) 


(ID 
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Figure  2.  Lossless  inductively  loaded  dipole. 

The  lossless  dipole  corresponds  to  a  collisionless  plasma. 


where 


Xe  “  “  4£2/Lu>2.  '12) 

If  we  have  a  random  distribution  of  dipoles,  we  must,  naturally,  con¬ 
sider  the  various  orientations  that  the  dipoles  can  assume.  As  shown 
in  figure  3  let  the  E  field  be  in  the  Z  direction.  For  every  random 
dipole  direction  0,  we  have  an  equally  probable  direction  -0.  We  pan 
therefore  pair  all  the  dipoles  oriented  in  the  +0  direction  with  dipoles 
•  oriented  in  the  -0  direction^  The  resultant  effective  dipole  moment  for 
a  dipole  is 

p(0)  =  -  2  cos  0  (£  21  cos  0)  E  e^ut/Lw2. 

The  average  value  for  a  random  distribution  of  directions,  pay,  is  then 

tt/2  /  tt/2 

Pav  ■  J  “4  cos20(£2E  e^a)t/Loj2)2Tr  sin  0  d©/ 2  f  2it  sin  0  d0 
0  7  o 

■  j  p  *  -  4£2E  e^Ut/3Lu2. 

Therefore,  we  have 


X  =  -  4&2/3Lw2. 
Aeav 
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Figure  5.  Two  dipoles  oriented  in  the  8  and  -0  directions. 
The  y  components  of  the  two  dipoles  cancel. 


Now  assuming  that  xmav  "  t'ie  index  of  refraction  of  the  artificial 
dielectric  would  be 

n  =  £  i  (1  -  4N£2/3e0L(D2,>  j  J1/2 

and  the  intrinsic  impedance  would  be 

n  =  ^yo/eoU  -  4N£2/3e0Lu2)]  l^2 

In  the  case  of  tenuous  plasma  (ionized  high-altitude  air)  we  have  an 
index  of  refraction  given  by 

n  =  (1  -  tfe2/eonui)2) */2  =  (1  -  3.18*  109i^/ u>2) 1  (14) 

where  N  is  the  electron  density  in  cm"3,  m  is  the  electron  mass,  and  e 
is  the  electron  charge.  Therefore,  if  we  want  to  simulate  a  tenuous 
plasma  corresponding  to  electron  density  iV,  we  can  determine  the  length 
21  and  the  inductance  L  of  our  inductively  loaded  dipoles  with  the 
relation 


Mi  =  1  Mi 

£qL  4  G0m 


2.39*10%  sec"2 
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N&2 

L 


2.11*  10“ 2iV  sec-2  F/tn 


(15) 


or 


where  N  is  the  density  of  artificial  dielectric  dipoles  in  meter-3,  £  is 
in  meters,  and  L  is  in  henrys. 

Next,  we  note  that  by  introducing  a  resistance  into  our  dipole  we 
can  simulate  the  imaginary  part  of  the  index  of  refraction  that  arises 
because  of  the  collisions,  v,  per  second  between  the  pl'-sma  electrons 
and  the  atmospheric  molecules.  Figure  1  shows  that  the  voltage  between 
the  dipole  spheres  is 

2tE  eJut  -  L  +  R  4^ 
dt*  dt 

where  as  before  we  assume  Lw  »  1/Cw;  Again,  we  can  assume  that  the 
charge  on  the  capacitance  between  the  two  spheres  varies  as 


Then 


2£E  e^Ut  =  LQ(-w2)e3Ut  +  hqjto  e^Wt 

and 

Q  =  -  2£E/(Lu>2  -  jRio) 
or 

q  =  -  2£E  ejut:/(Lu)2  -  jRw). 

The  dipole  moment  along  the  axis  of  the  inductive  dipole  in  figure  1  is 
then 

p  =  2£q  =  -  4£2E  e^Ut/(Lu2  -  jRw). 

We  see  by  analogy  with  the  purely  inductive  dipole 

Pav  =  -j  p  =  -  4£2E  e^Wt/3(Luj2  -  jRw). 
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Therefore,  we  have 


Xeav  =  "  4*2/30^2  - 


(16) 


and 


«  e0  [l  -  4£2N/3e0Lo32(l-  j R/toL)  j 


(17) 


■  e0 


1  - 


4  &2N 
3  eqL 

(f)2  +  *2 


+  j 


4  £2N  ,R  l 
3  cqL  L  oj 

d  2  - 

<|>  +w 


(18) 


Finally,  we  note  that  if  equations  (3)  and  (4)  are  satisfied,  we  have  a 
one  to  one  correspondence  between  equations  (1)  and  (19). 


2 . 2  Explicit  Consideration  of  Dipole  Capacitance 


Up  to  this  point  we  have  assumed  that  the  inductance  used  to  load  the 
short  dipole  dominated  the  behavior  of  the  dipole .  We  will  now  examine 
this  assumption  in  more  detail.  The  arguments  presented  here  are  basically 
due  to  Harrington.15  Consider  two  short  dipole  antennas  as  shown  in 
figure  4.  Antenna  1  is  excited  by  a  current  source  Ij;  antenna  2  is  loaded 
at  its  center  with  an  impedance  Zt,.  The  currents  and  voltages  of  the  two 
dipoles  can  be  related  by  the  standard  impedances  of  a  two-port  network: 


V1  “  Z11  *1  ~  z12  *2 
V2  =  Z21  *1  -  z22  I2* 


(19) 


The  voltage  V  can  be  expressed  in  terras  of  ZL  and  I2:  V2  =  ZL  I2.  Then 
we  can  express  the  current  passing  through  the  load  Zj,  in  terms  of  the 
current  source  driving  antenna  1;  the  transfer  impedance  Z2 i ;  the  driving 
point  impedance  of  antenna  Z,  Z22;  and  the  load  impedance,  Z^i 

Z21  I1 

lz  =  (zL  +  Z22T*  (20) 

We  next  make  use  of  the  reaction  concept,  first  introduced  by 
Rumsey.16917  Rumsey  difined  reaction  of  field  a  on  the  source  b,  {a,b}, 
with  the  following  integral: 
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Figure  4.  Two  short  dtlpole  antennas. 

I 

{a,b}  =  J'j'j  (v.3 -  Ha*Mb)  dv 

where  j}3  is  the  electric  current  density  of  source  b  and  M*5  is  the  mag¬ 
netic  current  density  of  source  b.  In  this  notation  the  reciprocity 
theorem  is: 

{a,b}  =  {b,a}. 

For  a  current  source  1^,  we  can  obtain 

fa.b}  =  f  Ea‘Ib  d£  =  Ib/  Ea*dfc  =  -  VaIb.  (21) 

Ine  transfer  impedance  in  equation  (19)  is  defined  by 


where  V^j  is  the  voltage  at  terminal  i  produced  by  the  current  source  at 
terminal  j.  In  equation  (21)  we  indicate  that 

{j,i}  =  -  Vyli, 
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J 


therefore,  we  obtain  by  eliminating  V^j  and  using  the  reciprocity  theorem 

z«  ’  -  Tfiffh  (22> 

We  next  remove  dipole  number  1  to  minus  infinity  so  that  the  electric 
field,  Ei,  at  dipole  number  2  is  a  plane  wave.  King18  shows  that  for  a 
short  dipole  of  length  B,  i.e.  for  (tt/X)  B  <  1/2,  we  can  assume  an 
induced  current  distribution  on  dipole  2  of  the  form: 

11  (z)  =  I2(l-fjz|).  (23) 

Equation  (22)  then  yields: 

ExB 

z12  =  2jY  =  Z21  (24) 

and  from  equation  (20)  we  obtain: 

Z21  1 1  E1B 

12  “  (2l  +  Z22)  "  2(ZL  +  Z22)* 

Equation  (23)  then  can  be  expressed  as 
.  ElB  , 

X2(z)  =  2(Zl  +  Z22)  (1  ‘  B  (25) 

The  conservation  of  charge  along  dipole  number  2  can  be  expressed 
as 


d  ijw 

£ - +  j  w  q(z)  -  0,  (26) 


where  q(z)  is  now  the  charge  per  unit  length,  and  therefore  the  charge 
j distribution  along  dipole  2  is: 


q(z) 


d  4^ 

dz 


(27) 


The  dipole  moment  along  antenna  2,  p,  is  given  by  the  following  integral: 


B 

2 

p  ~  f  q(a)  z  dz.  (28) 
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(29) 


Expressions  (25),  (26),  and  (28)  then  yield: 
-j  Ej  B2 

P  =  4w(j  wL  +  Z22) * 


If  we  allow  x  to  be  the  polarizability  of  the  dipole  we  have: 

P  c  E*x 
or 

_ -.1  B2 

X  4oj  (juL  +  Z22)  * 


(30) 


(31) 


Equation  (31)  is  essentially  the  same  as  the  result  we  obtained  with  the 
simple  quasi-static  model,  except  for  the  presence  of  the  driving  point 
impedance  Z22  in  the  denominator.  (The  factor  4  in  the  denominator  occurs 
because  we  have  not  capacitively  end-loaded  the  dipole  as  in  figure  1.) 

King19*20  derives  a  value  for  the  impedance  Z22  of  a  short  dipole: 


z22 


<f>2 

12 


j 


41n  (|)  -  6.78 


(32) 


Here  n  is  the  impedance  of  free  space,  120  ttO,  B  is  the  length  of  the 
dipole,  and  A  is  the  dipole  radius. 


Up  to  this  point  we  have  considered  the  incident  electric  field  E, 
to  be  parallel  to  the  short  inductively-loaded  dipole.  Of  course,  the 
current  and  charge  distribution  along  the  dipole  is  a  function  of  the 
orientation  of  the  receiving  antenna  with  respect  to  the  surface  of 
constant  phase  of  the  incident  electric  field.  King19  indicates  that 
for  a  short  receiving  antenna  (with  a  triangular  current  distribution) 
the  projection  of  the  electric  field  onto  the  antenna,  multiplied  by 
the  actual  half  length  of  the  antenna-  gives  the  emf  of  the  equivalent 
circuit.  That  is  equation  (24)  is  modified  to 

EXB 

z12  "  Z21  "  77o  cos  6  (33) 


where  0  is  the  angle  between  the  incident  electric  field  vector  and 
the  antenna.  For  example,  0=0®  when  the  antenna  lies  in  the  sur¬ 
face  of  constant  phase.  We  can  now  use  the  same  orientation  argu¬ 
ments  that  we  used  in  the  development  of  our  quasi-static  model  which 
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ignored  the  splf  impedance  Z22*  If  we  wish  to  calculate  the  polarization 
in  the  direction  of  the  incident  electric  field  when  the  Ej  vector  is  not 
parallel  to  the  dipole,  we  must  introduce  another  cos  0  and  integrate 
over  all  directions.  Therefore,  as  before 


P 


av 


and 
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av 


*av 


E 


or 


=  -  ,1  B2 

av  12a)(ju)L  +  Z02)  * 


(34) 


If  we  set  Z22  =  0,  B  =  2 Z,  and  end-load  the  dipole  so  that  the  current 
distribution  is  uniform  and  not  triangular,  we  obtain  the  same  result 
as  given  by  equation  (13). 18 


The  appearance  of  Z22  in  the  denominator  of  expression  (34)  is  in¬ 
consistent  with  our  desire  to  model  a  simple  Lorentzian  plasma.  Spe¬ 
cifically,  when  the  frequency  is  such  that 


wL 


41n(~)  -  6.78 


(?)  B 


(35) 


the  artificial  dielectric  will  pass  through  a  resonance.  We  can  only 
use  the  inductively-loaded  dipoles  to  simulate  a  plasma  if  10L  »  Z22» 
We  will  discuss  the  effect  of  Z22  in  a  subsequent  section.  An  insight 
into  the  limitations  imposed  by  the  dipole  self-impedance  can  be  ob¬ 
tained  by  calculating  xav  f°r  sorae  realistic  frequencies  and  dipole 
parameters. 


2 . 3  Mutually  Perpendicular  Dipole  Scatterers 

Our  arguments  presented  so  far  assumed  that  the  distance  between 
the  inductively-loaded  dipoles  is  great  enough  so  that  their  interaction 
can  be  neglected.  The  use  of  equation  (8)  is  based  on  the  assumption 
of  negligible  interaction  between  scattering  objects.  A  way  of  in¬ 
creasing  the  value  of  N  in  (8)  by  a  factor  of  3  and  still  not  violate 
the  assumption  of  negligible  interaction  of  neighboring  scatterers  is 
to  construct  each  scatterer  out  of  three  mutually-perpendicular 
inductively-loaded  dipoles.  By  using  essentially  quasi-static  argu¬ 
ments  it  is  easy  to  see  that  three  mutually-perpendicular  dipoles  with 
a  common  center  point  do  not  interact. 
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Consider  just  two  perpendicular  dipoles  as  shown  in  figure  5.  By 
symmetry  the  current  distribution  along  the  two  arms  of  dipole  1  is 
symmetrical  and  equal  in  magnitude.  Therefore,  the  charge  distribution 
along  the  opposite  arms  is  equal  in  magnitude  but  of  opposite  sign. 

Using  a  static  argument  the  electric  field  produced  by  the  charge  distri¬ 
bution  along  dipole  1  is  seen  to  be  perpendicular  to  the  direction  of 
dipole  2.  The  dipoles  therefore  do  not  interact.  Similar  reasoning 
applies  to  a  dipole  that  is  mutually  perpendicular  to  dipoles  1  and  2 
and  centered  at  the  intersection  of  dipoles  1  and  2.  These  quasi-static 
arguments  are  consistent  with  the  more  rigorous  arguments  presented  for 
skew-angled  dipoles  by  Richmond.21  If  we  substitute  three  mutually- 
perpendicular  dipoles  for  the  single  randomly-distributed  dipoles  we 
simply  multiply  the  average  value  of  the  polarizability  given  in  equation 
(34)  by  3  and  obtain 


Xav3 


_ zJL.Bg _ _ 

4to(j(oL  +  Z22)  * 


(36) 


By  constructing  our  scatterers  out  of  three  mutually-perpendicular 
inductively-loaded  dipoles  we  accomplish  more  than  simply  increasing  the 
number  of  scatterers  by  three:  the  polarizability  of  the  scattering 
centers  is  made  to  be  independent  of  scatterer  orientation.  To  be 
specific,  the  scatterer,  consisting  of  three  identical,  short,  mutually- 
perpendicular  dipoles  centered  at  the  origin  with  a  dipole  along  each 
axis,  can  be  described  with  the  following  simple  dyadic: 

A  A  A  A  A  A 

X  =  X  (i  i  +  j  j  +  k  k)  =  XI. 

The  dipole  moment  induced  by  an  arbitrary  field  E  is  then 
p  =  X  •  E 

=  X^  •  E  (i  sin  6  cos  <t>  +  j  sin  0  sin  <J>  +  k  cos  0)  (37) 

a  a 

=  X  E  (:'.  sin  0  cos  <5>  +  j  sin  6  sin  <i>  +  k  cos  0) 

but  (sin  0  cos  <t)2  +  (sin  0  sin  $)2  +  cos20  =  1.  Therefore,  the  magnitude 
of  the  dipole  moment  is  a  E  and  its  direction  is  in  the  direction  of  the 
electric  field  regardless  of  the  orientation  of  the  scatterer. 

In  order  to  gain  insight  into  the  foregoing  result  we  can  consider 
a  relatively  simple  situation  in  more  detail.  Figure  6  illustrates  two 
perpendicular  inductively-loaded  dipoles  in  the  plane  of  the  page.  Let 
the  E  vector  of  the  incident  plane  wave  also  be  in  the  plane  of  the  page. 
The  third  dipole  is  perpendicular  to  the  plane  of  the  page,  and  therefore 
is  perpendicular  to  the  E  vector  and  does  no£  interact  with  the  E  field. 

As  shown  in  the  figure,  the  Poynting  vector  S  of  the  incident  plane  wave 
is  also  in  the  plane  of  the  page.  The  voltage  developed  across  the  load 
of  dipole  number  1  is  then  proportional  to  E  cos  0  and  the  separation  of 
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■DIPOLE  ONE 


Figure  5.  Two  perpendicular  inductively-loaded  dipoles 
do  not  interact. 


charge  on  dipole  number  1  in  the  direction  of  the  incident  E  vector  is 
proportional  to  cos29,  that  is  P1,e||  =  K  cos20.  A  corresponding  argu¬ 
ment  (see  fig.  7)  leads  to  the  conclusion  that  the  dipole  moment  of 
dipole  number  2  in  the  direction  of  the  incident  field  E  is  given  by 
P2  e| I  =  K  sin20.  The  total  dipole  moment  in  the  direction  of  the  vector 
E  of  dipoles  1  and  2  is,  therefore,  given  by 


PK||  =  Pl,E||  +  P2,E|  |  =  K  COs20  +  K  Sitl20 

f38) 

.  K  =  . “I-5-25  - 

4u(ja)L  +  Z22)  ’ 


Let  us  now  consider  the  polarization  perpendicular  to  the  E  vector.  As 
shown  in  figure  8  the  voltage  developed  across  the  load  on  dipole  1  is 
proportional  to  E  cos  0  and  the  dipole  moment  perpendicular  to  the  E 
field  or  parallel  to  the  S  vector  ispijjj^-K  cos  0  sin  0.  As  shown 
in  figure  7,  since  the  voltage  across  tfie  load  on  dipole  2  is  proportional 
to  E  sin  0,  the  dipole  moment  component  of  dipole  two  perpendicular  to 
the  incident  field  is 
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Figure  6.  The  vector  is  the  dipole  moment  component  of 

dipole  1  that  is  parallel  to  the  incident  electric 
field. 


p0  „ .  =  K  sin  0  cos  0. 

Z  y  111 

So  for  the  total  polarization  perpendicular  to  1l  we  have 


Figure  7.  The  dipole  moment  of  dipole  2  can  be  resolved  into 

components  parallel  and  perpendicular  to  the  incident 
electric  field. 


In  summary,  we  see  that  the  rather  tedious  arguments  leading  to  equations 
(38)  and  (39)  are  consistent  with  the  result  of  equation  (37).  Of  course, 
as  the  length  of  the  dipoles  is  increased,  the  simple  trigonometric  de¬ 
pendence  for  the  current  or  charge  distribution  does  not  hold,  and  the 
polarization  is  not  independent  of  the  direction  of  the  incident  electric 
field. 


Coupling  of  Scatterers 


The  dipole  interaction  of  the  scattering  objects,  each  consisting  of 
three  mutually-perpendicular  and  centered,  inductively-loaded  dipoles, 
will  be  considered  next.  The  starting  point  for  the  discussion  of  the 
coupling  of  scattering  objects  in  both  real  and  artificial  dielectrics 
is  the  formula  of  Lorentz  for  the  field  within  a  hollow  spherical  cavity 
(centered  at  a  particular  molecule  or  scatterer)  cut  out  of  the  dielectric 
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Figure  8.  The  vector  p.  „  is  the  dipole  moment  component  of 

1  f&l 

dipole  1  that  is  perpendicular  to  the  incident  field. 
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(fig.  9).  The  displacement  vector  D  normal  to  the  dielectric  surface 

between  the  sphere  and  the  cavity  must  be  continuous;  therefore,  the 

magnitude  of  P  normal  to  the  cavity  surface  equals  the  negative  of  the 

induced  surface  charge  density  o±.  We  may  write  ci  =  -  jpj  cos  6.  i 

The  field  inside  the  cavity  due  to  C£  can  be  calculated  at  the  center  ( 

of  the  cavity  with  Coulomb's  law: 
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Therefore,  the  local  field  at  the  center  of  the  cavity  at  a  particular 
molecule  or  scatterer  is 

i 

E.  =  E+— -P  +  E2. 
loc  3eo 

Here  E  is  the  ordinary  macroscopic  field,  P  is  the  polarization  of  the 
dielectric,  and  E2  is  the  mean  field  due  to  the  matter  within  the  sphere 
itself  (exclusive  of  the  central  scatterer) .  The  magnitude  of  E2  will 
depend  on  the  distribution  of  scatterers  within  the  spherical  cavity. 
Lorentz  showed  that  if  the  molecules  are  arranged  in  a  cubical  array,  E2 
will  be  zero.  Lorentz  also  indicated  that  to  a  certain  degree  of  approx¬ 
imation  E2  -  0  for  such  isotropic  bodies  as  glass,  fluids,  and  gases. 

We  will  discuss  this  approximation  in  a  later  section. 

Assuming  that  E2  *  0,  we  obtain  the  following  value  for  the  local 
field  at  a  scatterers 

E1qc  =  E  +  P/3e0  (40) 

The  dipole  moment  of  a  scatterer  is  then  p  =  Xav3  ®loc»  ^ut> 

?  -  ?  N  ■  xav3  »  «loc 

where  N  is  the  number  density  of  scatters  per  unit  volume.  According 
to  the  usual  notation: 

e  E  =  e0  E  +  P. 

Therefore 

P  =  E  (e  -  eq) 
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and 
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E  ‘  *av3  E  Eloc  '  *av3  »  <E  +  E«c0) 
*)* 

E(ei-e0) 

■  Eav3S  <E+^— >' 


(41) 


This  is  the  Clausius-Mossotti  relation.  Equation  (41)  can  be  manipulated 
to  yield: 
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1  t  2  H  xaT3/3«, 

1  -  S  Xav3/3e„  ' 


(42) 


he  above  expression  is  only  valid  when  the  local  electric  field  E^oc  at  a 
scattering  obstacle  can  be  expressed  as 


hoc  =  E  +  P/(3e0)  (43) 


where  E  is  the  space  average  field  inside  the  artificial  dielectric. 


The  local  fields  for  simple  tetragonal  and  hexagonal  lattices  have 
also  been  calculated  by  a  number  of  workers.  Collin'*22  considers  a 
three-dimensional  array  of  y-directed  unit  dipoles  at  x  =  na,  y  =  mb, 
and  z  =  sc.  He  excludes  the  dipole  atx=y=z=0  and  calculates  the 
local  field  at  x  *  y  =  z  *  0.  The  potential  due  to  a  unit  dipole  at 

xq*  yo*  z0  ^ 
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4ir  gq  r3 


(44) 


where  r  =  [(x  -  x0)2  +  (y  -  yo)2  +  (Z  -  Zq)2]1/2.  For  the  array  of 
y-directed  dipoles  we  then  have  (fig.  10) 
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(45) 


where  the  prime  indicates  the  omission  of  the  n  =  m  *»  s  «  0  term.  The 
effective  polarizing  field  at  x  =  0,  y  =  0,  z  =  0  will  then  be  E  +  E^, 
where  E  is  the  external  applied  field  and  where 


i 


i 


i  . 

i  ' 


i 


i 


s  t* 


i 

i 


28 


CUBIC  LATTICE  a  =  b  =  c 
TETRAGONAL  LATTICE  a=c*b 
HEXAGONAL  LATTICE  a*b*c 

Figure  10.  Three-dimensional  array  of  indoctively-loaded 
dipoles  enclosed  in  styrofoam  spheres. 
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(46) 


Note  that  the  above  expression  is  not  only  evaluated  inside  the  Clausius- 
Hossotti  sphere,  but  over  all  space.  We  therefore  do  not  have  to  consider 
a  contribution  from  the  inside  surface  of  the  sphere. 


Collin  evaluates  expression  (46)  and  obtains 
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K0  (^>  +  K„ 


(47) 


where  Kg  is  the  so-called  modified  Bessel  function  of  order  zero.  Equa¬ 
tion  (47)  is  for  unit  dipoles.  In  our  case  of  inductively-loaded  dipoles 
we  would  have 
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or 
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The  above  equation  corresponds  to  a  cubic  lattice  if  a  =  b  =  c,  a  tetra¬ 
gonal  lattice  if  only  two  of  the  quantities  are  equal,  and  a  hexagonal 
lattice  if  none  of  the  quantities,  a,  b,  or  c  are  equal.  In  the  case  of 
the  cubic  lattice  we  find  that 
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which  is  to  be  compared  with  equation  (42) 
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2.5  Deviations  from  the  Clausius-Mossottl  Equation  Due  to  Randomness 
of  Scatterers 

J.K.  Kirkwood2  3  >2t*  has  calculated  the  deviations  that  one  might  expect 
from  equation  (40)  in  the  case  of  a  random  distribution  of  molecules  con¬ 
sisting  of  hard  spheres.  Specifically,  Kirkwood  considered  a  molecular 
interaction  potential  of  the  form 

W(r)  =»  0  _<  r  <_  a 

(52) 

=  0  r  >  a 

where  a  is  the  diameter  of  a  molecule. 

Kirkwood’s  model  is  not  wholly  applicable  to  an  artificial  dielectric 
formed  by  dumping  spherical  scatterers  into  a  container  because  the 
scatterers  do  not  necessarily  assume  a  completely  random  distribution. 

Eor  example,  depending,  among  other  factors,  on  how  smooth  they  are,  the 
spheres  can  arrange  themselves  into  various  specific  lattices  or  combina¬ 
tions  of  lattices.  The  different  lattices  can  even  correspond  to  dif¬ 
ferent  dipole  number  densities. 

It  is  interesting  to  note  that  statistical  mechanical  averaging 
based  on  the  molecular  interaction  equation  (52)  yield:  a  value  of  the 
permittivity  that  is  not  a  function  of  the  thermodynamic  temperature  T. 

This  follows  because  the  value  of  the  exponential  term  exp  {-W(r)/kT]  is 
either  0  or  1.  In  a  gas  or  fluid  the  averaging  can  be  over  time;  in  our 
static  collection  of  scattering  dipoles  the  averaging  would  have  to  be 
over  an  ensemble  of  different  containers  each  filled  with  a  aggregation 
of  spherical  dipole  scatterers.  Kirkwood's  approach  can  be  summarized 
as  follows:  As  in  the  lattice  case,  the  local  field  is  obtained  by 
summing  over  the  potentials  due  to  individual  dipoles.  We  can  say  that 
the  effective  polarizing  field  E^oci  at  scatterer  i  will  be  given  by  E  +  Ej, 
where  E  is  the  external  applied  field  and  E^  is  given  by  the  sum  over  N 
individual  dipoles: 


-y.  -y. 

E.  =  -  E  T.,  .p, 

1  .  .  =dk  rk 

k=l 

where  the  dipole-interaction  dyadic  is  given  by 

rp  =  _i _ !_  T  ,  rik  rik 

=ik  4tts0  3  =  2 

rik  rik 


(55) 


The  polarization  of  an  individual  scatterer  is  given  by 
=  e°='^loci  =  e0“*(E  +  E±) 

“  e0oE  -  e0a  E 

We  therefore  have  N  simultaneous  equations 


P±  +  ae0  E  Tik-Pk  =  ae0E  i  =  1 ,...N 

that  would  have  to  be  solved  in  order  to  obtain  Eioci  for  each  of  the 
individual  scatterers  in  the  distribution  of  N  scatterers  in  volume  V. 
We  can  average  the  N  equations  and  obtain 


p  +  ae0  E  =  ae0£ 

Kirkwood23  >2lt  incrtv.jced  the  following  fluctuation  term 


N 


5  '  (Stk'?k  -  Sik-fk5 

i^k 


and  obtained 


P  +  ae0  E  T±k*P  +  n<xeQ  =  oe02 


If  we  can  assume  that 


ll2*Pl  *  il2*Pl 

the  fluctuation  n  is  equal  to  zero  and 
p  +  ae0  p* E  =  ae0£. 

Converting  to  the  equivalent  scalar  equation  we  obtain: 
aeoE 


P  = 


1  +  aeoU  Tik) 


(56) 


(57) 


(58) 


(59) 


(60) 


(61) 
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and 


D  =  eE  =  e0E  +  P 


(62) 


.  N  - 
eoE  +  ~  P 


=  £qE  + 


eO  v  “E 


1  +  ae0  (E  Tik) 


(63) 


The  dipole  interaction  term  can  be  averaged  over  the  volume  v  to  yield: 
1  N 


E  T 


ik 


also 


3  v  ’ 


(64) 


a  = 


_av3 

£0 


(65) 


and  therefore 


or 


=  e0/l  + 


N 
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(66) 


(67) 


Equation  (67)  is  to  be  compared  with  equations  (51)  and  (42) .  After 
considerable  mathematical  manipulation  and  a  number  of  approximations, 
Kirkwood  managed  to  evaluate  equation  (57)  for  a  substance  consisting 
of  a  random  distribution  of  hard  noninteracting  spheres  of  diameter  a: 


,1_  _5,  N2a2D 

Sb  "  48v;  v 


(68) 
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where 


t.  ..  P-TiNa3  _  ,  ,ira3N 
b  =  -y  -  4iV  (— ) 


=  4  (volume  of  spherical  scatterers) . 

To  obtain  a  rough  estimate  of  the  validity  of  using  the  Clausius-Mossotti 
equation  for  a  random  distribution  of  styrofoam  ball  dipole  scatcerers, 
we  can  compare  the  value  of  n  given  by  equation  (68)  with  the  value  of  23 


,  .  1  Na 
1_aETik=  1  +  3  v 


For  the  experimental  situation  discussed  in  Section  2: 


n 


1  -  a  £  T 


ik 


.22 


and  we  might  expect  the  Clausius-Mossotti  equation  to  be  a  reasonable 
first  approximation.  Actually,  as  we  have  pointed  out,  Kirkwood's 
statistical  model  is  not  necessarily  applicable  to  our  aggregation  of 
dipole  scatterers  because  the  dipole  scatterers  in  the  artificial  die¬ 
lectric  do  not  necessarily  assume  a  random  distribution.  The  validity 
of  the  Clausius-Mossotti  equation  as  applied  to  our  artificial  dielectric 
should  really  be  tested  empirically. 

3.  COMPARISON  OF  EXPERIMENT  WITH  THEORY 


In  order  to  orient  ourselves  in  regard  to  the  parameters  such  as  the 
dipole  load  Impedance  and  Z22  we  will  use  equations  (36)  and  (42)  to 
design  an  artificial  dielectric  with  some  actual  experiments  in  mind. 

The  frequency  range  considered  will  be  dictated  by  the  band  pass  char¬ 
acteristics  of  the  type  2300  waveguide.  The  type  2300  waveguide,  the 
largest  standard  waveguide,  has  internal  dimensions  of  23.0  x  11.5  in., 
a  lower  cutoff  frequency  of  2.56  x  108  H,  and  a  recommended  upper  fre¬ 
quency  limit  of  4.9  x  108  H.  The  frequencies  in  this  range  are  much 
higher  than  the  values  usually  quoted  for  nuclear  EMP  pulses,  but  we  will 
consider  designing  an  artificial  dielectric  to  operate  at  these  high  fre¬ 
quencies  because  of  the  convenience  of  using  waveguide  techniques  in 
measuring  the  properties  of  an  artificial  dielectric. 

Let  us  consider  the  following  series  of  parameters  as  an  orientation 
exercise.  We  begin  with  an  artificial  dielectric  constructed  out  of 
scatterers  consisting  of  three  mutually-perpendicular  inductively-loaded 
7-in.  long  dipoles  loaded  with  a  number  of  different  values  of  inductance. 
Let  the  dipole  stems  have  a  radius  of  0.1  cm.  We  choose  the  density  of 
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scatterers  to  be  264  per  cubic  meter.  This  density  corresponds  to  an 
aggregation  of  48  scatterers  distributed  throughout  a  rectangular  con¬ 
tainer  with  dimensions  23  x  11.5  x  41.5  in.  Equations  (36)  and  (42) 
can  be  employed  to  calculate  the  dielectric  constant  of  the  artificial 
dielectric  as  a  function  of  frequency  in  the  band-pass  region  of  the 
type  2300  waveguide.  Figure  11  shows  values  of  the  real  part  of  the 
permittivity  calculated  with  equations  (36)  and  (42)  for  a  frequency 
range  between  300  and  500  MH.  The  four  curves  correspond  to  dipole 
load  inductances  of  0.4,  0.6,  0.9,  and  1.4  yH.  The  magnitude  of  the 
radiation  resistance  of  the  dipoles  is  very  small  compared  to  the  values 
of  the  magnitude  of  the  inductive  loads;  therefore,  the  imaginary  com¬ 
ponents  of  the  permittivity  are  very  small. 


If  we  want  to  model  a  Lorentz  plasma  in  which  the  electrons  have  a 
collision  frequency  v  t  0,  the  load  impedance  of  our  dipoles  must  have  a 
resistive  component.  Equation  (36)  then  becomes: 


-  i  B2 

*av3  4u(jwL  +  R  +  Z22) 


(69) 


where  R  is  the  value  of  the  resistance.  In  figures  12a  to  12d  we  present 
values  of  the  real  part  of  the  permittivity  calculated  with  equations  (69) 
and  (42)  with  L  =  G.4  pH.  The  series  of  four  figures  shows  the  effect  of 
different  values  of  the  resistive  component  R  on  the  magnitude  of  the 
resonance  occurring  when  juL  =  Z22*  figure  12a,  corresponding  to  a  purely 
inductive  load,  emphasizes  the  resonance  behavior  of  the  artificial  die¬ 
lectric  produced  when  the  load  inductance  L  is  in  resonance  with  the 
capacitive  impedance  of  the  dipole.  As  can  be  expected  the  Q  of  the 
resonance  is  reduced  by  including  a  resistance  in  the  dipole  load  im¬ 
pedance. 


Figure  12a,  The  real  component  of  relative  permittivity  calculated 
with  equations  (69)  and  (42).  The  value  used  for  L  is 
.4  pH  and  the  value  used  for  R  is  On.  The  other 
parameters  are  the  same  as  those  employed  to  obtain 
figure  11. 


Hr  x  10' 


Figure  12b.  The  real  component  of  relative  permittivity  calculates 
with  equations  (69)  and  (42).  The  value  used  for  L  is 
.4  pH  and  the  value  used  for  R  is  50ft.  The  other 
parameters  are  the  same  as  those  employed  to  obtain 
figure  11. 


L  =  .4  pH 
R  -  100  n 


Hr  x  10“ 

Figure  12c.  The  real  component  of  relative  permittivity  calculated 
with  equations  (69)  and  (42) .  The  value  used  for  L  is 
.4  pH  and  the  value  used  for  R  is  100ft.  The  other 
parameters  are  the  same  as  those  employed  to  obtain 
figure  11. 


Figure  12d.  The  real  component  of  relative  permittivity  calculated 
with  equations  (69)  and  (42) ,  The  value  used  for  L  is 
.4  yH  and  the  value  used  for  R  is  200ft.  The  other 
parameters  are  the  same  as  those  employed  to  obtain 
figure  11. 


In  our  first  experiment  we  searched  for  the  resonance  predicted  by 
equation  (36)  and  shown  in  figure  12a.  The  - cperimental  technique25 
was  straightforward  and  is  depicted  in  figure  13.  The  ratio  of  the 
incident  electric  field  Ej  to  the  reflected  electric  field  ER  was  ob¬ 
tained  by  measuring  the  standing  wave  ratio  with  a  slotted  line.  The 
artificial  dielectric  sample  shown  in  figure  13  consists  of  an  aggrega¬ 
tion  of  48  of  our  7  in.  styrofoam  balls  each  containing  3  mutualiy- 
perpendicular  dipoles.  The  balls  were  enclosed  in  a  rectangular  box 
measuring  23  x  11.5  x  41.5  in.  The  results  of  the  experiment  are  shown 
in  figure  14.  It  is  comforting  to  note  that  there  is  a  large  reflection 
at  approximately  the  frequency  of  the  resonance  shown  in  figure  12a. 

In  a  second  series  of  experiments  we  measured  the  permittivity  of 
our  artificial  dielectric  utilizing  the  well-known  shorted  waveguide 
technique  described  in  great  detail  by  von  Hippel26.  Very  briefly,  the 
technique  consists  of  measuring  the  null  point  of  a  standing  wave  in  a 
shorted  waveguide.  A  sample  of  length  l  is  then  inserted  into  the 
shorted  waveguide  and  the  shift  in  the  standing  wave  null  is  noted.  The 
index  of  refraction  of  the  dielectric  sample  of  length  £  ce:i  then  be 
calculated  in  terms  of  %  and  the  shift  in  the  standing  wave  null  point. 
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Figure  13.  Experimental  setup  used  in  determining  the  resonance 
frequency  of  the  artificial  dielectric.  The  sample 
consisted  of  48  scatterers  or  "molecules."  The 
V.S.W.R.  or  voltage  standing-wave  ratio  was  determined 
with  a  slotted  line. 


FREQUENCY  IN  MECAHERTZ 


Figure  14.  Experimental  standing-wave  ratio  corresponding 
to  the  setup  shown  in  figure  13. 


Our  experimental  values  of  the  index  of  refraction  and  some  calculated 
values  are  shown  in  figure  15.  The  measurements  were  of  two  different 
artificial  dielectric  sample  sizes;  one  sample  consisted  of  an  aggrega¬ 
tion  of  24  scatterers  measuring  23  x  11.5  x  20.75  in.  in  volume.  The 
length  of  the  sample  was,  £  =  20.75  in.  The  other  sample  consisted  of 
an  aggregation  of  48  scatterers  in  a  rectangular  volume  measuring 
13  x  11.5  x  41.5  in  In  this  case,  £  =  41.5  in. 

The  difference  between  the  theoretical  index  of  refraction  curve 
corresponding  to  L  =  0.4  pH  and  R  =  0  and  the  measured  index  of  re¬ 
fraction  is  about  30  percent.  A  possible  source  of  the  disparity  be¬ 
tween  experiment  and  theory  is  that  we  used  the  dimension  of  the  card¬ 
board  box  that  enclosed  the  scatterers  to  obtain  £.  The  effective 
length  of  the  very  granular  artificial  dielectric  sample  is  most 
probably  not  the  size  of  the  box  containit  the  scatterers.  The 
boundary  of  our  artificial  dielectric  is  somewhat  nebulous  and  needs 
more  study. 

Further  experimental  measurements  of  the  electromagnetic  character¬ 
istics  of  our  inductively-loaded  dipoles  ere  now  under  consideration. 
These  new  measurements  would  involve  the  irradiation  of  a  monopole  that 
is  immersed  in  an  expanse  of  artificial  dielectric  medium.  The  new 
experimental  approach  would  overcome  the  limitations  imposed  by  the 
relatively  small  size  and  narrow  frequency  band  pass  of  a  waveguide. 

The  measurements  could  also  be  carried  out  at  lower  frequencies  which 
are  closer  to  the  actual  nuclear  EMP  spectrum.  As  the  frequency  is 
lowered,  the  size  of  the  scatterer  compared  to  the  wavelength  would 
decrease.  That  is  to  say,  the  artificial  dielectric  would  have  less 
granularity  and  the  boundary  of  the  artificial  dielectric  would  become 
more  well  defined. 

If  we  construct  a  cylindrical  monopole  of  height  H-  it 'would  have 
a  resonance  or  maximum  skin  current  at  frequency  f  given  by 


If  the  monopole  were  then  submerged  into  an  artificial  dielectric  with 
an  index  of  refraction  n,  the  resonance  of  the  monopolr  would  change  to 


The  influence  of  the  artificial  dielectric  of  the  resonance  character¬ 
istics  and  skin  currents  of  other  shapes  such  as  a  sphere  could  also  be 
studied. 


j 


40 


4.  SUMMARY  AND  CONCLUSIONS 


We  first  introduced  a  simple  quasi-static  model  of  the  artificial 
dielectric  consisting  of  inductively-loaded  dipoles  which  (1)  ignored 
the  dipole  capacitance,  and  (2)  assumed  that  the  dipole-dipole  inter¬ 
action  .  ”jld  be  neglected.  These  two  simplifying  assumptions  yielded 
an  expression  for  the  permittivity  of  the  artificial  dielectric  that 
behaved  in  a  manner  surprisingly  analogous  to  the  permittivity  of  a 
simple  Lorentzian  plasma.  It  was  noted  that  the  granularity  of  fluc¬ 
tuations  of  the  artificial  dielectric  could  be  decreased  by  constructing 
the  dipole  scatterers  out  of  three  mutually-perpendicular  scatterers. 

Then  the  effect  of  the  capacitance  of  the  inductively-loaded  dipole 
scatterers  was  examined  and  limits  set  on  the  assumption  that  the  in¬ 
ductance  dominated  the  behavior  of  the  inductively-loaded  dipole.  Next, 
the  effect  of  dipole-dipole  interactions  on  the  artificial  dielectric 
was  considered.  A  relatively  rigorous  expression  for  a  cubic,  tetra¬ 
gonal  or  hexagonal  lattice  was  obtained.  An  estimate  of  the  effect  of 
a  random  distribution  of  scatterers  on  the  value  of  the  permittivity 
predicted  by  equation  (42)  was  also  obtained. 

It  is  interesting  that  the  more  exact  expressions  for  the  dielectric 
constant  of  the  artificial  dielectric,  which  considered  both  the  effect 
of  the  capacitive  dipole  impedance  and  the  dipole-dipole  interaction, 
did  not  model  the  behavior  of  a  Lorentzian  plasma  as  well  as  the  original 
simplified  approach.  Nevertheless,  the  artificial  dielectric  did  have 
an  index  of  refraction  less  than  one  over  a  relatively  broad  spectrum. 

We  are  now  looking  into  the  possibility  of  improving  correspondence 
between  the  frequency  dependence  of  our  artificial  dielectric  and  a 
Lorentzian  plasma.  The  approach  being  investigated  is  based  on  the 
substitution  of  a  more  complicated  one-terminal  pair  in  place  of  the 
simple  inductive  dipole  load  used  in  our  present  scatterers.  By  syn¬ 
thesizing  a  suitable  one-terminal  pair  out  of  a  number  of  lumped  im¬ 
pedances  we  hope  to  improve  the  correspondence  between  the  frequency 
dependence  of  the  artificial  dielectric  and  a  Lorentz  plasma. 

The  experimental  examination  of  the  theoretical  expressions  for  the 
dielectric  constant  of  our  artificial  dielectric  was  constrained  by  the 
dimensions  and  frequency  characteristics  of  available  waveguide  equip¬ 
ment.  The  predictions  of  equations  (36)  and  (42)  for  the  dielectric 
constant  of  the  artificial  dielectric  were  found  to  be  in  fair  agree¬ 
ment  with  experimental  results. 
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